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Laser-plasma accelerators have become compact sources of ultrashort electron bunches at energies 
up to the gigaelectronvolt range thanks to the remarkable progresses made over the past decade. A 
direct application of these electrons bunches is the production of short pulse x-ray radiation sources. 
In this letter, we study a fully optically driven x-ray source based on the combination of a laser- 
plasma accelerator and a plasma wave undulator. The longitudinal electric field of a laser-generated 
plasma wave is used to wiggle electrons transversally. The period of this plasma undulator being 
equals to the plasma wavelength, tunable photon energies in the ten keV range can be achieved with 
electron energies in the 100-200 MeV range. Considering a tens terawatt class femtosecond laser 
system, undulators with a strength parameter K ~ 0.5 and with about ten periods can be combined 
with a laser-plasma accelerator, resulting in several 10"'^ emitted x-ray photons per electron. 



Laser wake field acceleration (LWFA) has been one of 
the most developed topic of high intensity laser-plasma 
interaction during the past ten years. It has lead to num- 
ber of theoretical and experimental efforts resulting in 
remarkable results. In LWFA, electrons are accelerated 
from rest to relativistic energies thanks to strong longi- 
tudinal electric field (lOO's GeV/m) of a large amplitude 
plasma wave driven in the wake of an intense femtosec- 
ond laser pulse. Electron beams currently produced are 
quasi-monoenergetic [IHSj with an energy spread down 
to below the percent H], an energy up to the GeV [S], a 
charge in the 10-100 picocoulomb range, a divergence of 
a few mrad, a duration in the few femtoseconds range [6] , 
and a normalized emittance in the tt. mm. mrad range. 

In the interaction regimes of laser wakefield electron 
acceleration, there usually occurs the emission of in- 
tense and collimated femtosecond x-ray bursts. Sev- 
eral mechanisms, based on the radiation from relativis- 
tic electrons from plasmas, have been studied as sources 
of femtosecond x-rays: the non-linear Thomson scatter- 
ing of a laser beam by plasma electrons [71 [5] , the Be- 
tatron oscillations of accelerated electrons in the bub- 
ble regime [9,, the Thomson backscattering of a laser 
beam by laser-accelerated electrons [10], the oscillations 
of laser-accelerated electrons in a periodical structure of 
magnets (the conventional meter-scale undulators used 
in synchrotron facilities) [HI [12] . These high-brightness 
compact and ultrashort x-ray sources can potentially 
have many applications in fundamental research as well 
as in the industrial or medical domains. However, even 
if these sources can already cover a spectral range from 
a few eV (conventional undulators) to a few 100 keV 
(Thomson backscattering) there exists a gap in the ten 
keV range where none of these sources is efficient. Only 
the Ka source emits in this energy range but it is a fully 
isotropic radiation and has a duration of a few hundreds 
femtosecond. The production of a radiation beam with 
a few femtosecond duration at this energy is therefore 
particularly important for applications. Indeed, most of 
the applications require a source with an angstrom wave- 



length and a femtosecond duration, respectively the typ- 
ical interatomic distances and time scales of fundamental 
processes in matter. 

In this letter, we discuss the possibility to produce 
a fully optically driven femtosecond and tunable x-ray 
source in the ten keV range by using a laser-generated 
plasma wave as an undulator for a laser-accelerated elec- 
tron bunch. The concept of the plasma wave undulator 
was originally proposed in the context of realization of 
compact short wavelength free electron lasers from con- 
ventional accelerators [TS] [14]. Indeed, such an undu- 
lator has a very small period compared to conventional 
undulators T5; and allows to produce short wavelength 
free electron radiation with moderate electron energies. 
Further studies on the plasma wave undulator were per- 
formed a few years later |16H19j . discussing in more de- 
tails the classical spontaneous emission. Electron trajec- 
tories as well as the incoherent radiation they emit were 
analytically or numerically calculated using ideal linear 
plasma waves. In the following, we will study the utility 
of the plasma wave undulator for laser-plasma acceler- 
ators. The interest of this study is motivated by the 
active research on femtosecond x-ray sources and by the 
fact that laser-plasma accelerators have become efficient 
and reliable at the 100 MeV energy level. We will show 
that, at these electron energies, high brightness and tun- 
able undulator type radiation can be produced in the 
relevant ten keV range. Our description does not rely on 
ideal linear plasma wave. Instead, we first propose a con- 
cise analytical analysis of the nonlinear one-dimensional 
plasma wave undulator for estimation of the radiation 
properties, and then fully describe the realistic plasma 
wave generation using quasi-static particle code simula- 
tions . These realistic plasma waves are then used to 
calculate numerically electron trajectories and the emit- 
ted x-ray radiation. The consideration of realistic plasma 
waves is particularly important since for the considered 
amplitudes they become non linear, with possibly damp- 
ing and breaking of the plasma wave, and strong mod- 
ification of the driving laser pulse along its propagation 
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FIG. 1. Principle of the plasma undulator source combined 
with a laser-plasma accelerator. An electron bunch is created 
and accelerated at the interaction of an intense femtosecond 
laser with a millimeter scale plasma. Once accelerated elec- 
trons propagate and are wiggled in a plasma wave. Their 
relativistic oscillation motion results in the emission of an x- 
ray beam. 



in the plasma. All these effects are taken into account in 
the used quasi-static particle code WAKE [20*. 

A plasma wave corresponds to the coupled propagation 
of the electron density rie and the scalar/electric potential 
V. For non-linear multi-dimensional plasma waves, the 
charge current je and the magnetic field B also appear in 
their description. Several routes exist to produce plasma 
waves. They have been explored in order to produce in- 
tense electric fields for electron acceleration in plasmas. 
Large amplitude plasma waves can be generated through 
the beating of two copropagating long laser pulses satis- 
fying the condition uj2 — ui = ojp [wp = (nee^/meo)^^^ 
is the plasma frequency, uti and uj2 the laser frequencies 
of each laser pulse] , the self- modulation instability [a sin- 
gle long laser pulse breaks in a pulse train of duration 
T ~ Ap/c, where Ap = 27rc/wp is the plasma wavelength] 
and the resonant case where the laser pulse has a duration 
T ~ Ap/c. These methods have been demonstrated to be 
efficient for the generation of electric fields E = — 
in the lOO's GeV/m range (it has to be compared with 
the highest achievable electric field by conventional ap- 
proaches, which is in the lO's MeV/m range). In the fol- 
lowing, we will consider the generation of plasma waves 
by a resonant excitation. 

The schematic of the x-ray source studied is presented 
m Fig. [l] Elect rons are first accelerated to relativistic 
energies in a laser-plasma accelerator (along the x direc- 
tion) and then propagate in a plasma wave traveling in 
the z direction. As they propagate in the plasma wave, 
they experience successively an electric field directed to- 
ward the z positive values and toward the z negative 
values. As a consequence they oscillate in the z direction 
with a significant transverse acceleration. Since both the 
plasma wave and the electron bunch travel at almost the 
speed of light, the period of the electron motion is equal 
to the plasma wave period, i.e. approximately the plasma 



wavelength (which is 10.6 /im for an electron density of 
lO^^cm"'^). The result of this relativistic oscillatory mo- 
tion is the emission of an x-ray beam by the electron 
bunch. Experimentally this can be realized by using two 
synchronized laser pulses. The first laser pulse, directed 
along the x direction, is used to create the plasma ac- 
celerator, producing a relativistic electron bunch. The 
second laser pulse, propagating in the transverse direc- 
tion z, used to create the plasma wave undulator. 

The radiation produced in the plasma undulator is a 
relativistic moving charge radiation and its features de- 
pend on the electron motion. In order to derive esti- 
mations for the electron orbits driven by a plasma wave 
undulator, we consider a one dimensional plasma wave 
described by the normalized electric potential (/>(^) = 
eV/mc^ only depending on the comoving variable ^ — 
z — Vpt (quasi-static approximation) , Vp being the plasma 
wave phase velocity and supposedly equal to the laser 
group velocity. 

We assume an electron from a laser-plasma accelerator 
with an initial energy £ = jiUic^ (7^ is the initial rela- 
tivistic factor) and a velocity directed along the x axis. 
The electron enters transversally in the one-dimensional 
(ID) plasma wave potential 0(0. The equation of the 
test electron motion is: 



dp 



(1) 



where all quantities are normalized by the choice m ^ c = 
e = Wp = 1. The Hamiltonian describing the dynamic of 
the test electron in the plasma wave is: 



(2) 



where P — p — a is the canonical momentum. In our 
case the canonical and kinematic momentum p are equal 
because the vector potential a = Q. T-L does not depend 
on the x and y coordinates. This implies the conservation 
of px and Py. Furthermore it depends on t and z only 
through ^ = z — j3pt :^ z — t. This gives the following 
constant of motion: 



^ — 4> — Pz — const. 



(3) 



Using the initial conditions for the test electron: x = 

^/^l — 1, = 0, the above 
1 + p^, we 



y = z = Py = Pz Px = ^Ti (\> 

constants of motion and the relation 7^ 
obtain: 



27i li 



(4) 



where we develop expressions using the approximation 
^ 7j which is always satisfied (the normalized poten- 
tial cannot exceed a few units in realistic plasma waves, 
whereas 7^ is superior to 100 for our case). 
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The electron motion consists of a drift in the x direc- 
tion combined with an oscillation in the transverse direc- 
tion z. The plasma wave acts as an undulator with 
a period A,; ~ Xp. This period depends on the plasma 
density and can be controlled. In practical units, it is 
given by Xu[^J.ra] = 3.34 x 10^''' 



'ne[cm"3]. A second 
relevant parameter characterizing the plasma undulator 
is its strength parameter K. It is defined as the prod- 
uct of the electron energy 7 by the maximum angle of 
the trajectory ip (with respect to the x axis): K = 'yijj. 
Because V = dz/dx^^^ ~ dz/dt^^^ = _p2max/7, from 
Eq. Q the strength parameter K of the plasma wave 
undulator reads 



K = -fil: =p^ 



(5) 



where 4>o is the amplitude of the plasma wave poten- 
tial (/)(^). In the linear limit Oq ^ 1, the plasma wave 
reads — •\/7raQ(wpr/4) exp(— a;pr2/4) sin(fcp^), for a 
gaussian laser pulse with duration r at l/e^ in intensity, 
and if we approach the non-linear regime, we can reach 
K = ^o<l. 

The radiation emitted by an electron traveling in the 
plasma wave undulator is a moving charge emission de- 
scribed by the general expression giving the radiation 
emitted by an electron, in a direction of observation n, 
as a function of its position, velocity and acceleration 
along the trajectory [21j : 
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dt 



This expression represents the energy radiated within 
a spectral band duj centered on the frequency w and a 
solid angle dJl centered on the direction of observation 
n. Here, f{t) is the electron position at time t, /3 is 
the velocity of the electron normalized to the speed of 
light c, and /? — d/3/dt is the usual acceleration divided 
by c. The analysis of this general expression results in 
analytical expressions describing the radiation features 
as a function of the parameters 7^, A„ and K. In the 
laboratory frame, the radiation emitted by a relativistic 
electron at a given time is directed in the direction of its 
velocity v, within a cone of typical angle I/7. The angu- 
lar distribution of the radiation depends on the electron 
orbit. For a trajectory confined in the plane {x, z), which 
is the case for a ID plasma wave, the radiation is emit- 
ted in a cone of half angle 9 = 1/^ ior K < 1, while 
for i<r > 1 the angular distribution is increased in the 
direction of the oscillation z with a half angle of if/7. 
The electron motion being periodic with a spatial period 
A„, the spectrum of the radiation consists on the fun- 
damental wavelength A = ^{1 + /2 + j^9'^) and its 
harmonics, where 9 is the observation angle with respect 
to X. For the undulator limit K <^ 1, i.e small transverse 



amplitudes of oscillation, the spectrum contains only the 
fundamental wavelength, and harmonics start to appear 
when K 1. In the wiggler limit K ^ 1, the spec- 
trum contains many closely spaced harmonics extending 
up to a critical wavelength Ac = ^^^2 ■ For the case of 
the plasma wave K < 1 and the radiation is produced 
in the undulator regime at the fundamental wavelength 
and possibly a few harmonics. The spectral thickness of 
the fundamental wavelength and its harmonics depend 
on the number of oscillation periods N. It is given by 
AA„/A„ — l/{nN) where the index n indicates the har- 
monic order. Note that the spectrum is broadened with 
a decreasing energy when observed off axis. The total 
number of emitted photons depends on the electron en- 
ergy 7, the K parameter and the number of oscillations 
N . Over one oscillation period and for K <^ 1, it is given 
by Nx = ^aK^ (at the mean energy {huj) = j'^hc/Xu) 
PTj where a = e^/AiTeofic is the fine structure constant. 
Finally, the duration of the emitted radiation pulse is 
sensibly equal to the duration of the electron bunch. In 
the case of laser-plasma acceleration, the duration of the 
electron bunch can be as small as a few femtoseconds [6] . 

For example, if we consider a 150 MeV electron travel- 
ing in a plasma wave undulator with K — 0.5 and a pe- 
riod A„ = 10 /im the produced radiation has an on-axis 
fundamental energy Ex ~ 21 keV. The half angle of the 
x-ray beam divergence is ~ 3 mrad and the expected 
number of photons produced per electron and per period 
2 is ~ 4 X 10~^. Assuming 10^ electrons and 10 periods the 
total number of photons is ~ 4 x 10^. This photon en- 
ergy range is not covered by the existing laser produced 
x-ray sources. In addition with the present methods of 
Wakefield acceleration, the electron energy can typically 
be tuned from 50 to 200 MeV [12]) which corresponds to 
x-ray energy from ^ 2 to 37 keV. 

The above description assumes an idealized situation 
where the plasma wave is ID and the electron motion 
is strictly periodic. A more precise description based on 
numerical simulations is presented in the following. 

The generation of the plasma wave by a resonant ex- 
citation has been calculated using the Wake code [20j . 
Wake is a two dimensional, moving window particle code 
describing laser-plasma interaction in the short pulse rel- 
ativistic regime. It allows to calculate the propagation 
of a relativistic laser pulse in an underdense plasma and 
its wakefields. To speed up the calculation time com- 
pared to PIC (particle in cell) simulations. Wake assumes 
that the laser pulse does not evolve significantly during 
the electron transit time through the laser pulse (quasi- 
static approximation), and it uses averaged equation of 
motions for the particles (average over the laser period) 
so that it is only necessary to resolve the plasma time 
scale, and not the laser time scale. As a consequence, 
the code can not describe high energy electrons moving 
in the forward direction with 7 > uj^/wp, where loq is the 
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FIG. 2. The map of the field Ez obtained fi"om the Wake code 
is displayed in (a). It corresponds to the wakefield after 0.5 
mm of propagation of the laser pulse in the gas jet of density 
rie = lO'^^'cm"^. The profile -Ez(C) at y = is shown in (b). 



laser frequency. 

The wakefield maps obtained from Wake are then used 
to calculate the orbits of test electrons in the plasma wave 
undulator. From the electron trajectories, we compute 
the x-ray radiation properties using Eq. Q. 

We have considered a laser pulse, propagating in the 
z direction, focused on a gas jet with a density of Ue — 
lO^^cm"-^. The focal spot is chosen to be elliptical, with 
a waist larger in the x direction than in the y direction in 
order to increase the overlap of electrons with the plasma 
wave. The laser waist equals 11 /xm in the y direction, 
the pulse duration is 30 fs at 1/e^ in intensity, and the 
amplitude of the normalized vector potential of the laser 
pulse is ao — 1.2. After 0.5 mm of propagation in the 
plasma, the Wake simulation [run in the two dimensions 
{z,y)] shows that the amplitude of the plasma wave at- 
tains (/)o — 0.5 , as displayed on Fig. [2] The maps of the 
wake fields E^, Ey and B^, after 0.5mm of propagation in 
the plasma, characterize the plasma wave undulator. To 
take into account the finite dimension of the laser pulse 
and of the wake fields in the x direction, we apply a gaus- 
sian profile with a waist of 44 /im to the wake fields for 
this direction. 

We have injected an electron beam of energy £ — 
150 MeV and of transverse size 7 ^m, propagating in 
the X direction, in the plasma wave undulator charac- 
terized by the fields obtained from Wake, and we have 
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FIG. 3. x-ray angular profile of the plasma wave undulator 
radiation, in J/sr and per electron. 




FIG. 4. The on-axis spectrum {9z — 6y = 0) of the plasma 
wave undulator radiation is displayed in (a), in J/eV/sr and 
per electron. The integrated spectrum over angles is shown 
in (b), in J/eV and per electron. 



calculated the produced radiation using Eq. (|6|. The 
electrons emit a narrow x-ray beam directed in the x 
direction, with a divergence of a few mrad. The angu- 
lar profile of the radiation is shown on Fig. [3] Figure 
|4] displays the spectra of the radiation emitted on axis 
as well as integrated over the angular distribution. It is 
in agreement with the properties of the undulator radia- 
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tion, with an angle-dependent wavelength and a narrow 
bandwidth at a given angle of observation. The central 
photon energy and the full width at half maximum band- 
width of the on-axis spectrum are respectively 17 keV 
and 22%. The total number of emitted photons is 10~^ 
per electron. The peak brightness associated to the on- 
axis spectrum of Fig. |4] is on the order of 10^° — 10^^ 
photons/0. l%bandwidth/s/mm^/mrad^, considering 10^ 
electrons, a source size of 5-10 microns and a duration 
of a few femtoseconds. These results show good agree- 
ment with the estimations obtained with the analytical 
ID study. 

In conclusion, the combination of a laser-plasma ac- 
celerator with a plasma undulator offers the possibility 
to produce a compact tunable source of ultrashort x- 
ray radiation in the ten keV range. This source would 
be tunable in a wide spectral range thanks to both the 
tunability of the undulator period with the plasma den- 
sity and the electron energy. This x-ray source exhibits 
the properties of undulator radiation: narrow bandwidth 
spectrum and high brightness. 

The plasma wave undulator approach presents similar- 
ities with the Thomson scattering, and their advantages 
and disadvantages will be briefly discussed. In the Thom- 
son scattering geometry a laser pulse scatters off the elec- 
tron beam from the laser-plasma accelerator. Electrons 
oscillate in the laser pulse and emit radiation that can 
rapidly attain very high photon energies thanks to two 
successive Doppler shifts. The oscillation period is given 
by A„ ~ Ao/(l — /3cos0), where Aq is the laser wavelength 
and 9 the angle between the laser and the electron beam 
directions (co-propagating for 6 = 0), and the strength 
parameter of the oscillation is K = gq where ao is the 
normalized amplitude of the laser pulse vector potential. 
If we consider the highest quality electron beams from a 
laser-plasma accelerator produced to date, their energy is 
in the hundred MeV range. At these energies, the Thom- 
son scattering is more appropriate to produce energetic 
photons, in the 100 keV or MeV range. Indeed, the pro- 
duction of x-ray radiation at 10 keV using electrons at 
200 MeV requires a co-propagating geometry where the 
laser and the electron beam make about a 10 degrees an- 
gle. This complicates the experimental setup due to the 
positioning of the focusing optics. In addition, the laser 
scattering off the electron beam has to travel into the few 
millimeter gas jet, can diffract and lose energy. Both the 
Thomson scattering and plasma wave undulator sources 
can have a strength parameter K of the same order of 
magnitude, < 1, but it is easier to reach a higher number 
of oscillation periods with the Thomson scattering, and 
so a higher number of photons and brightness. 

The plasma wave undulator appears to be an interest- 
ing scheme because its period in the ten microns range 
is adapted to the electron energies which are currently 
produced. It allows to reach the 10 keV range using 



electrons at a few hundreds MeV. In addition, the un- 
dulator period and therefore the energy of the radiation 
can be easily tuned by varying the plasma density, keep- 
ing the electron beam constant. Even if the two schemes 
can therefore be used to access the ten keV range and 
produce femtosecond x-ray radiation, the Thomson scat- 
tering looks more appropriate to access higher energies, 
the hundred keV range, while the plasma wave undulator 
is suitable to access the ten keV range. 

In this letter we have chosen to discuss the properties 
of this source as it could be produced with 10-100 TW 
short pulse laser systems. The numerical simulations, in 
agreement with analytical predictions, show that a radi- 
ation beam, coUimated within a few mrad, at energies 
in the ten keV range and with ^ 10~^ photons / elec- 
tron can be produced. Assuming electron bunches with 
60 pC charge, the source can deliver ^ 10^ photons / 
pulse. The source size is equal to the dimension of the 
electron bunch within the undulator. This depends on 
the distance between the electron source and the undu- 
lator. Typically, a source size on the order of 5 /im can 
be expected. The flux of this source is not as high as the 
betatron source or the source, however this source 
combines the advantages of narrow band spectrum and 
tunability, the coUimation of a few mrad and the ultra- 
short pulse duration of a few femtoseconds. With more 
laser energy, a major improvement would be to increase 
the number of periods of the plasma wave undulator, by 
focusing on a larger waist in the direction of propagation 
of electrons a;, and by generating a plasma wave consist- 
ing of many plasma periods. Assuming a plasma wave 
undulator with 25 periods and with K ~ 0.5 the source 
could reach the level of ^ 10^^ photons / electron with 
an on-axis spectral bandwidth of a few percents. The use 
of future high repetition rate laser systems would allow 
to improve the averaged brightness of this source. 

We acknowledge Professor P. Mora for providing the 
Wake code. 
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